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1.0  INTRODUCTION 


I! 


Digital  techniques  have  opened  new  alleys  in  the 
field  of  ionospheric  sounding.  System  reliability  and  data 
accuracy  have  improved;  new  parameters  are  measured:  ampli- 
tude, phase,  incidence  angle,  polarization  and  Doppler  shifts; 
data  storage  on  digital  magnetic  tape  makes  computer  processing 
easy . 

Air  Force  Geophysics  Laboratory's  Ionospheric  Dynamics 
Branch  operates  two  Digisondes,  one  in  Goose  Bay.  Labrador  and 
the  second  on-board  a KC-135  aircraft.  We  have  modified  the 
aircraft  Digisonde  by  introducing  the  technique  of  spectral 
integration.  Phase  coherent  integration  is  a powerful  tech- 
nique that  significantly  increases  the  system  sensitivity  as 
compared  to  power  integration.  The  new  Processing  Controller 
that  was  integrated  in  the  aircraft  Digisonde  makes  a Discrete 
Fourier  Transform  (DFT)  in  each  frequency-range  bin  thus 
carrying  out  phase  coherent  integration  for  each  Doppler 
shifted  spectral  line.  The  Doppler  shifts,  of  course,  are 
introduced  by  the  moving  ionosphere  and/or  the  moving 
aircraft . 

To  make  the  Goose  Bay  observations  available  to  the 
scientific  community,  we  have  scaled  the  hourly  vertical  iono- 
grams  according  to  URSI  specifications.  The  quarter-hourly 
ionograms,  both  vertical  and  oblique  backscatter,  are  avail- 
able on  magnetic  tape  and  paper  printout. 


2 . 0  DIGISONDE  SYSTEMS 
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2.1  Goose  Bay  System 

■Because  of  the  large  dynamic  variations  in  signal 

field  strengths  a digital  automatic  gain  control  was  incor- 
porated in  the  Digisonde  in  addition  to  the  already  existing 
' time  dependent  gain  control.  The  latter  changes  the  receiver 

i gain  in  steps  of  10  dB  in  the  evening  and  in  the  morning  to 

compensate  for  the  diurnal  variation  in  D region  absorption. 
The  new  automatic  gain  control  takes  for  each  transmitted  fre- 
; quency  the  range  bin  with  the  largest  amplitude  as  decision 

maker.  If  in  the  initial  phase  of  the  integration  time  the 
amplitude  is  smaller  than  90%  of  the  end  value  the  receiver 
gain  is  changed  by  one  or  two  10  dB  steps.  The  modified  re- 
ceiver gain  is  stored  in  the  preface  (position  23)  of  the 
output  data. 

The  two  Granger  wide  band  RF  amplifiers  were  made 
compatible  by  introducing  new  wide-band  transformers.  One  of 
the  plate  circuit  delay  lines  was  replaced.  Tests  and  cali- 
brations were  performed  to  demonstrate  proper  operation  of  the 
amplifiers. 

2.2  Aircraft  System 

The  phase  coherent  integration  in  the  original  Digi- 
sonde could  not  be  applied  if  a significant  Doppler  shift  was 

imposed  on  the  signal.  The  Doppler  frequency  is  given  by 
I 

1 = ' 

where  k is  the  wave  vector,  specifying  the  looking  angle,  ^ 
is  the  aircraft  velocity  and  is  the  velocity  vector  of  the 
ionosphere . 


The  new  Processing  Controller  that  is  now  incorpor- 
ated in  the  aircraft  Digisonde  is  not  merely  overcoming  this 
problem  by  using  spectral  integration,  but  by  measuring  the 
Doppler  frequencies  it  makes  available  a new  parameter  in 
ionosonde  observations.  The  features  of  the  Processing  Con- 
troller are  described  in  Appendix  A. 


2.3  Data  Replay  System 

The  Digicoder  of  AFGL's  data  replay  system  was  modi- 
fied so  as  to  allow  for  greater  formatting  flexibility.  The 
printout  can  use  three  different  fonts  now:  the  original  Digi- 
coder font  with  the  6 by  4 matrix,  a double  sized  Digicoder 
font  with  a 12  by  8 font  and  64  level  BCD  font  with  a 12  by  8 
font  (Figure  1). 
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3.0  DATA  PROCESSING 
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3.1  lonogram  Evaluation 

To  obtain  a good  survey  of  the  diurnal  and  seasonal 
variations  of  the  arctic  ionosphere  at  Goose  Bay,  the  hourly 
vertical  ionograms  from  March  1975  to  February  1976  were 
scaled  according  to  URSI  definition.  In  addition,  we  scaled 
the  top  frequencies  and  virtual  heights  of  oblique  Es  and  F 
echoes,  generally  assumed  to  indicate  the  location  of  the 
northern  edge  of  the  ionospheric  trough.  The  following  23 
parameters  were  scaled:  foE,  foFl,  foF2,  foFOB*,  foEs,  foEs2*, 
frEs,  fbEs,  fminFOB*,  fol,  foT,  h'E,  h'Es,  h'Es2*,  h'F,  h'F2, 
h'FOB*,  kpFl,  MUF(3000)F2,  MUF(3000)F1,  Type  Es  and  Type  Es2 
(*  = oblique  echoes). 

In  Appendix  B we  show  the  diurnal  variations  of  the 
monthly  median  values  of  h'E,  h'Es,  h'F,  h'F2  and  foE,  foEs, 
foFl,  foF2.  As  an  indicator  of  the  sporadic  E activity  Appen- 
dix B also  presents  the  foEs  distribution  curves.  It  should 
be  noted  that  the  night  E,  a relatively  thick  layer  at  E region 
altitude,  that  frequently  forms  in  the  arctic  ionosphere  due  to 
particle  ionization,  has  been  included  in  these  foEs  distribu- 
tion curves. 

The  individual  hourly  values  are  kept  on  file  and 
are  available  on  request. 

3.2  Seminar  in  lonogram  Reduction 

A comprehensive  course  in  scaling  of  digital  iono- 
grams was  taught  at  the  University  of  Lowell  to  personnel  from 
AFGL.  The  course  used  as  text  the  "URSI  Handbook  of  lonogram 
Interpretation"  by  Piggot  and  Rawer.  As  a complementary  aid 
we  had  prepared  the  "Succinct  Guide  to  the  Reduction  of  Digi- 
tal Ionograms"  which  is  tailored  to  the  Digisonde  ionograms 


and  emphasizes  the  special  conditions  occurring  in  the  arctic 
ionosphere.  This  little  guide,  attached  as  Appendix  C to  this 
report,  is  very  helpful  for  new  personnel  trained  on  ionogram 
scaling,  although  it  is  not  in  any  way  a complete  ionogram 
scaling  manual. 
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ABSTRACl  I 

I 

1 

The  digital  ionospheric  sounding  system  "Digisonde 
12 BPS"  explores  ionospheric  structure  and  dynamics  by  exploiting 
all  observables  of  the  reflected  electromagnetic  wave:  range, 
amplitude,  phase,  Doppler,  incidence  angle  and  wave  polarization. 

The  Digisonde  operates  in  two  complementary  modes  of  observations: 
the  lonogram  Mode  with  full  range  and  frequency  display  but  lim- 
ited resolution  in  Doppler  and  incidence  angle,  and  the  Doppler- 
Drift  Mode  operating  with  a limited  number  of  frequency-range  bins 
but  full  resolution  in  Doppler  and  incidence  angle.  This  tech- 
nique reduces  the  volume  of  magnetically  recorded  output  data  to 
a manageable  size.  Special  cigital-analog  hybrid  techniques  are 
used  for  the  presentation  of  the  multi-parameter  ionograms  and 
the  sky  maps. 

j 

I 

INTRODUCTION 

Remote  sensing  of  the  atmosphere  is  a continuous  re- 
quirement for  the  understanding  of  our  environment  and  for  the 
warning  of  natural  and  man-made  events  that  affect  man's  life. 
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The  Wide  Band  Satellite  program  conducted  by  the  U.S.  Defense 
Nuclear  Agency  for  the  investigation  of  the  aurora  and  equatorial 
ionosphere  required  ground  based  equipment  for  continuous  moni- 
toring of  the  ionosphere  with  regard  to  vertical  and  horizontal 
electron  density  distributions  and  to  the  motions  of  ionospheric 
irregularities.  The  explanation  of  scintillation  and  spread  F 
phenomena  are  prime  goals  of  this  program.  The  tasks  of  an  iono- 
sonde  are  to  study  and  monitor  the  ionosphere  by  radio  means. 

The  research  objectives  of  the  ionosonde  demand  special  capabili- 
ties and  adequate  flexibility  from  the  instrument,  requirements 
which,  in  general,  are  expensive.  In  contrast,  the  most  impor- 
tant requirement  for  a monitoring  ionosonde  is  low  cost,  since 
only  a dense  network  of  sounders  can  produce  a global  picture  of 
the  ionosphere.  It  appears,  therefore,  that  two  sets  of  iono- 
sondes  are  needed,  one  for  research  and  one  for  monitoring  pur- 
poses. Such  a diversification,  however,  would  be  very  unsatis- 
factory unless  the  "monitoring  ionosonde"  (M.I.)  can  easily  be 
expanded  to  a "research  ionosonde"  (R.I.),  and,  also  important, 
the  R.I.  can  fulfill  the  operational  requirements  of  an  M.I. 

Our  answer  to  the  existing  need  was  the  development  of 
the  Digisonde  (Bibl  et  al,  1970;  Phillips,  1974)  an  advanced  dig- 
ital ionosonde.  Since  1969  Digisondes  were  used  for  both  iono- 
spheric monitoring  (AFCRL's  "Geophysics  and  Space  Data  Bulletin"; 
lonospharen-Institut  Breisach  "lonospharen-Daten" ; Institut  Royal 
Meteorologique  de  Belgique  "Observations  lonospheriques " ) and  for 
research  (Bibl  et  al , 1975;  Reini.ch,  1973;  Philbrick  et  al,  1973 
and  1974;  MacLeod  ct  al,  1973;  Reinisch  and  MacLeod,  1973).  The 
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initial  difficulty  of  the  digital  sounder  was  the  high  cost  of 
the  components  and  of  the  peripheral  recording  equipment . In  the 
last  three  years  very  versatile  digital  integrated  circuits  came 
on  the  market  at  low  prices , especially  the  so-called  read-only- 
memories  (ROM's)  and  the  microprocessors.  This  allowed  us  to 
build  the  new  Digisonde  128PS  at  moderate  costs  with  all  the 
capabilities  of  an  R.I. 

OBSERVABLES  FOR  lONOSPHERICALLY  REFLECTED  RADIO  SIGNALS 


The  ideal  ionosonde  must  be  able  to  measure  all  param- 
eters that  characterize  the  reflected  radio  wave.  Some  new  iono- 
sondes  were  recently  described  in  the  literature  emphasizing 
either  one  wave  parameter,  for  instance  the  phase  (Hammer  and 
Bourne,  1976)  for  accurate  virtual  height  determination  or,  low 
costs  (Morgan  and  Pratt,  1975).  It  appears  that  only  digital 
techniques  can  succeed  in  determining  all  the  wave  parameters. 

Fortunately,  the  list  of  parameters  is  not  too  long: 
group  travel  time  or  h'  (1),  amplitude  (2),  phase  (3),  precise 
frequency,  i.e.  Doppler  offset  (4)  from  transmitted  frequency, 
incidence  angle  (5),  and  wave  polarization  (6),  i.e.  ordinary  or 
extraordinary  echo.  The  curvature  (7)  of  the  wave  front  may  be 
of  interest  for  certain  experiments.  To  obtain  a complete  pic- 
ture of  the  structure  and  of  the  motions  in  the  ionosphere,  it  is 
desirable  to  measure  parameters  2 through  7 as  functions  of  the 
virtual  height  h'  and  the  transmitted  frequency  f over  large 
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ranges  of  h'  and  f,  and  with  good  resolution,  i.e.  small  incre- 
ments in  range  and  frequency.  But  the  amount  of  data  would  be- 
come excessive  unless  digital  preprocessing  and  multiplexed  out-  * 

put  formatting  reduce  the  volume  of  output  data  to  manageable 
size.  The  argument  could  be  made,  that  rather  than  outputting 
the  measured  data  one  could  immediately  extract  the  physical  in- 
formation which  is  of  interest  in  the  very  moment , and  not  store 
the  measured  data.  This  approach  may  lead  to  erroneous  results 
without  the  possibility  of  data  reverification  and  reinterpreta- 
tion. The  solution  to  the  apparent  dilemma  is  to  preprocess  the 
data,  reduce  the  redundancies,  and  flag  each  data  point  with  a 
parameter  identification.  The  appropriate  identification  tech- 
niques will  be  explained  in  the  next  sections. 

The  simultaneous  monitoring  of  the  entire  five-dimen- 
sional complex  amplitude  space  (range,  frequency,  Doppler,  azi- 
muth and  zenith  angle)  is  neither  feasible  nor  advisable.  In- 
stead, the  Digisonde  128PS  operates  in  two  complementary  modes, 
the  ionogram  and  the  drift  mode.  By  complementary  we  mean  that  a 
combination  of  the  two  data  sets  generated  in  the  two  modes,  con- 
tains all  relevant  information  of  the  five-dimensional  space.  In 
the  ionogram  mode,  128  or  256  range  bins  (virtual  heights)  are  sampled 
in  a selected  band  of  transmitted  frequencies,  with  full  ampli- 
tude resolution  (1/2  or  1 dB),  but  limited  information  about  the 
incidence  angles  and  Doppler  frequencies.  The  wave  polarization 
is  identified  for  all  frequency-range  bins  (FRB).  In  the  drift 
mode,  only  a limited  number  of  I’KB's  are  selected  and,  of  course. 
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these  will  be  FRB's  with  echo  signals.  A priori  assumptions 
about  the  pulse  shape  of  the  echoes  and  about  the  selection  of 
specific  frequencies  during  the  ionogram  (Wright,  1975)  are 
avoided  in  order  not  to  bias  the  two  basic  data  sets.  Table  1 
lists  the  wave  parameters  that  the  Digisonde  observes  and  the 
techniques  that  are  used  to  measure  these  parameters. 

SYSTEM  CONFIGURATION 


The  Digisonde  is  housed  in  two  standard  19"  racks  (Fig- 
ure 1),  one  holding  the  actual  ionosonde,  the  other  the  10  kW 
wide-band  pulse  transmitter.  Figure  2 shows  the  block  diagram  of 
the  sounding  system.  Besides  the  10  kW  TRANSMITTER  there  are 
only  four  chassis:  the  PROCESSING  CONTROLLER,  the  TRANSCEIVER, 
the  ANTENNA  SWITCH,  and  the  MICROCOMPUTER.  The  digital  output 
data  are  recorded  on  one  7 -track  magnetic  tape  drive,  for  both 
ionogram  and  drift  measurements.  The  ionograms  are  displayed  in 
real  time  on  a high  intensity  plasma  display  and  simultaneously 
recorded  on  an  electrostatic  printer.  A digital-analog  hybrid 
method  is  applied  for  the  data  presentation  (Patenaude  et  al, 
1973),  as  illustrated  in  Figure  3.  The  small  optically  weighted 
number  font  gives  the  data  printout  the  appearance  of  an  intensity 
modulated  analog  picture  while  it  retains  the  digital  information. 
In  Figure  3,  the  amplitude  values  are  given  in  4 dB  increments. 
After  ten  data  lines  a preface  line  is  inserted  containing  date, 
time  (hour,  minutes,  seconds)  and  all  program  parameters  (pulse 
repetition  rate,  receiver  gain,  etc.). 

u 5 
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The  PROCESSING  CONTROLLER  is  the  brain  of  the  Digisonde 


system.  All  control  functions  branch  out  from  here  to  the  rest 
of  the  system,  and  the  entire  digital  data  preprocessing  is  per- 
formed here,  including  the  discrete  Fourier  transform.  A small 
keyboard,  mounted  on  the  front  panel,  enables  selection  of  the 
operational  parameters.  Data  may  also  be  entered  on  a teletype 
terminal  which  is  preferable  when  a typed  record  of  the  entry  is 
desired . 

The  TRANSCEIVER  is  the  analog  signal  processor  which  is 
controlled  by  digital  signals  from  the  PROCESSING  CONTROLLER. 

The  digitally  synthesized  transmitter  frequency  is  fed  through  a 
two  stage  active  filter  with  400  kHz  bandwidth  that  is  kept  con- 
stant over  the  entire  frequency  range  from  .5  to  40  MHz.  The 
same  bandpass  filters  are  alternately  used  as  tuned  input  stages 
for  signal  reception,  protecting  the  first  mixer  from  high  level 
f off-frequency  interferers  that  would  otherwise  cause  undesirable 

cross-modulations,  and  as  tuned  output  stages.  The  1.  local 
oscillator  is  kept  2 MHz  off  the  transmitter  frequency  resulting 
in  a 2 MHz  first  intermediate  frequency  (IF)  which  is  then  con- 
verted to  the  final  IF  of  255  kHz.  The  IF  bandpass  filter  is 
composed  of  ten  uncoupled  LC  circuits  with  individual  bandwidths 
of  70  kHz,  resulting  in  an  overall  bandwidth  of  15  kHz. 

The  real-time  post-processing  is  performed  in  the 


MICROCOMPUTER.  Its  main  task  is  the  calculation  of  the  sky  maps 
which  give  an  immediate  picture  of  the  horizontal  structure  of 
the  ionosphere.  Cross-correlation  of  the  complex  spectra  obtained 


for  the  different  antennas  determines  the  incidence  angle  for 
each  spectral  component.  A true  height  algorithm  converts  the 
virtual  heights  of  the  echo  traces  in  the  ionogram  into  a vertical 
electron  density  profile  which  is  plotted  onto  the  ionogram  on 
the  Plasma  Display.  The  MICROCOMPUTER  can  also  be  used  off-line 
for  the  processing  of  ionogram  and  drift  tapes. 


MULTI-PARAMETER  lONOGRAMS 


Operational  Programs.  Three  different  ionogram  pro- 
grams are  stored  in  the  Digisonde,  two  of  which  can  be  programmed 
by  the  keyboard  or  the  computer,  and  the  third  is  maintained  in  a 
non-volatile  FROM  ( f ield -programmable  read -only-memory ) . The 
ionogram  rate  can  be  selected  out  of  a set  of  12  ionogram  se- 
quences, extending  from  1 to  120  ionograms  per  hour  (Table  2). 
lonograms  can,  of  course,  also  be  started  by  manual  or  computer 
commands.  The  most  important  programmable  operational  parameters 
are  listed  in  Table  3.  Most  entries  in  the  Table  are  self- 
explanatory,  but  some  items  need  to  be  discussed.  Item  #3,  for 
instance,  lists  either  equidistant  frequency  increments  or  a dual 
step  size,  for  example  10  and  90  kHz.  The  dual  mode,  selected  for 
phase  measurements,  increases  the  group  height  resolution: 
h*  = (c/4ti)  A((i/Af  , 

(c  = speed  of  light,  A^/Af  = change  of  phase  with  frequency). 

The  large  frequency  steps  provide  for  high  resolution,  and  the 
small  steps  resolve  the  2n  ambiguity.  For  the  10  and  90  kHz 


mode,  the  resolution  in  virtual  height  is 

fih'  = (c/4tt)  6(A^)/Af  = (c/4ir)  2ir/16 • 9 •lO*  = 100  m 
and  the  2v  ambiguity  range  is 

Ah'  = (c/4ir)  2^/10"  Hz  = 15  km, 

which  is  easily  resolved  by  the  amplitude-range  measurement.  If, 
and  only  if  any  of  the  dual  frequency  increment  modes  are  se- 
lected, will  the  phase  be  recorded. 

The  range  increments  (item  #4)  can  be  varied  over  a 
large  set  of  possibilities  from  0.5  to  10.5  km.  The  Digisonde 
samples  two  times  128  ranges  and  the  two  range  groups  can  form 
one  large  range  window,  or  two  non-abutting  windows.  To  obtain  a 
better  group  height  resolution  in  the  E-region  than  in  the  F- 
region  a bi-linear  scale  can  be  selected.  For  instance,  the  E- 
region  can  be  sampled  with  1.5  km  increments  and  the  F-region 
with  3.0  km.  The  phase  coding  of  the  transmitted  radio  waves 
makes  it  possible  to  identify  echoes  in  a later  pulse  period. 

The  phase  code  sequence  is  a quasi-random  series  of  180®  phase 
shifts  applied  to  the  HF  of  consecutive  pulses.  This  "keys"  the 
individual  transmitter  pulses  and  permits  the  use  of  high  pulse 
repetition  rates.  The  coding  also  helps  to  suppress  quasi- 
coherent  interferers  in  the  course  of  the  spectral  integration 
process . 

Doppler  and  Incidence  Angle.  Spectral  analysis  is  per- 
formed separately  for  each  FRB.  Since  it  would  be  impossible  to 
record  a complete  spectrum  of  say  128  lines  for  each  bin,  the 
spectrum  is  limited  to  16  lines,  arranged  symmetrically  around 
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zero  Doppler  shift.  The  Doppler  frequencies  together  with  the 
incidence  angle  measurements  help  to  interpret  the  ionograms  in 
tei’ms  of  the  structure  and  the  motions  in  the  ionosphere,  and  we 
expect  that  our  measurements  will  provide  new  inputs  to  the  phe- 
nomenon of  spread  F.  Phase  coherent  integration  without  spectral 
resolution,  i.e.  assuming  zero  Doppler  frequency,  would  render 
the  system  insensitive  to  moving  reflectors  like  the  auroral 
ionosphere  or  chemical  releases.  Without  phase  coherent  integra- 
tion, on  the  other  hand,  the  ionosonde  would  be  useless  in  all 
densely  populated  areas  with  high  prevailing  radio  interference 
levels.  An  escape  to  ever  higher  transmitter  power  is  not  the 
answer  in  such  circumstances,  since  that  would  turn  the  ionosonde 
into  an  unacceptable  source  of  interference  for  other  radio 
services . 

The  azimuthal  component  of  the  incidence  angle  is  mea- 
sured for  each  FRB  with  a resolution  of  30®.  For  this  measure- 
ment, ten  of  the  antennas  of  the  receiving  array  are  summed  via 
switched  delay  lines.  The  delays  are  determined  in  such  a way 
that  the  zenith  angle  is  kept  equal  to  15®  while  the  azimuth 
angle  is  scanned.  Different  zenith  angles  can  be  chosen  by  re- 
placing the  set  of  delay  lines.  For  each  of  the  selected  azimuth 
angles  the  signals  are  phase  coherently  integrated  for  each  FRB, 
respectively.  These  measurements  are  virtually  simultaneous 
since  the  typical  scanning  time  is  30  msec.  After  all  the  sam- 
ples of  a given  transmitter  frequency  are  taken,  the  azimuthal 
incidence  in  each  height  bin  is  determined  by  a differential 


I 

technique.  The  six  amplitudes  (in  case  of  60®  resolution),  ob- 
tained for  the  six  azimuthal  directions  are  compared  and  the  max- 
imum amplitude  together  with  the  corresponding  azimuthal  identi- 
fier are  retained  in  the  output.  Other  combinations  of  incidence 
angles  can  be  selected,  as  vertical-north-north/ east , or  verti- 
cal-south-south/west, etc.  The  differential  method  is  very  sen- 
sitive even  when  the  beam  width  formed  by  the  antenna  array  is 
not  small.  As  a matter  of  fact,  the  beam  width  should  be  big 

enough  to  avoid  deadspots  on  the  sky.  For  this  reason  we  se- 

lected the  inner  ten  antennas  of  the  receiving  array  for  the 
ionogram  observation. 

The  same  differential  technique  is  applied  for  all  the 
other  ionogram  parameters.  After  the  set  of  spectral  lines  is 
calculated  the  maximum  line  together  with  the  line  number  is  de- 
termined for  output  storage.  Each  of  the  ten  inner  antennas  sep- 
arates the  ordinary  and  the  extraordinary  modes  by  using  crossed- 
loop  antennas.  At  latitudes  for  which  the  quasi-longitudinal  ap- 
proximation holds,  the  two  loops  are  added  and  subtracted  via  a 

in/ 2 phase  shifter;  at  the  magnetic  equator  the  two  perpendicular 
loops  are  used  independently.  Again,  for  each  FRB  both  antenna 
outputs  are  compared  (after  integration);  the  larger  amplitude  is 
stored  and  decides  whether  the  FRB  receives  an  0 or  an  X flag. 


COMBINATION  OF  lONOGRAM  PARAMETERS 


AND  THE  lONOGRAM  OUTPUT  FORMAT 

Not  all  the  wave  parameters  are  recorded  with  maximum 
resolution  in  each  ionogram.  The  combination  of  parameters  de- 
pends on  the  type  of  observational  program  that  is  carried  out. 

On  the  other  hand  it  would  be  unwise  to  change  the  data  output 
format  each  time  the  measuring  program  is  changed.  Based  on  the 
experience  we  had  during  the  past  ten  years  with  the  older  Digi- 
sonde  systems  we  opted  for  a standard  format  for  the  recording  on  ' 

magnetic  tape.  The  format  is  indeed  the  same  as  the  format  of 
the  earlier  Digisondes.  For  each  transmitted  frequency  128  nine- 
bit  data  points  are  recorded  together  with  an  identifying  pref- 
ace. The  nine  data  bits  comprise  six  amplitude  bits  and  three 
status  bits  that  describe  the  selected  wave  parameters.  In  some 
cases  the  amplitude  resolution  is  decreased  to  five  bits  in  order 
to  increase  the  status  word  to  four  bits.  An  unambiguous  identi- 
fication of  the  amplitude  and  status  bits  is  always  contained  in 
the  preface.  Seven-track  digital  tape  drives  were  chosen  for  re- 
cording, since  only  a minority  of  computer  facilities  can  handle 
nine-track  tapes. 

It  was  stated  earlier  that  the  Digisonde  collects 
either  128  or  256  range  samples.  In  the  latter  case  the  total 
range  is  divided  into  two  range  windows.  Before  storing  the  data 
in  the  128  output  array,  the  differential  method  is  applied 
again,  selecting  for  each  FRB  the  larger  amplitude  from  the  two 
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range  windows;  the  respective  window  number  is  stored  in  the 
status  word.  The  probability  of  two  echoes  residing  in  the  same 
bin  number  in  windows  1 and  2 is  generally  less  than  5%,  unless 
severe  range  spread  conditions  prevail,  in  which  case  the  opera- 
tor can  reduce  the  repetition  rate. 

The  maximum  number  of  data  integration  channels  for 
each  FRB  is  24,  limited  by  the  memory  capacity  in  the  PROCESSING 
CONTROLLER  and  the  speed  of  the  microcircuits.  Hence  not  all  the 
parameters  can  simultaneously  be  measured  with  maximum  resolu- 
tion. A total  of  32  parameter  combinations  (Table  4)  can  be  se- 
lected by  varying  the  number  of  spectral  lines  (1  to  16),  range 
windows  (1  or  2),  polarizations  (1  or  2)  and  incidence  angles 
(1  to  6).  The  pulse  repetition  rate  and  the  number  of  time  sam- 
ples used  in  the  calculation  of  the  Fourier  series  determine  the 
spectral  frequencies.  The  operator  must  merely  select  which  fre- 
quencies he  wants  to  retain  for  recording  by  using  a table.  This 
is  a decision  between  good  Doppler  resolution  or  wide  Doppler 
range.  One  must  recognize  that  the  scimpling  period  doubles  when 
two  antenna  configurations  are  sampled  in  alternation,  and  in- 
creases proportionately  for  more  configurations.  In  turn,  the 
frequencies  of  the  selected  spectral  lines  reduce 
proportionally . 

lonogram  Display.  Real-time  display  and  printout  of 
the  digital  ionograms  is  a necessity  for  an  R.I.  and  at  least 
desirable  for  an  M.I.  A special  technique,  originally  developed 
for  the  Digisonde,  uses  optically  weighted  fonts  for  the  printout 
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of  the  ionograms  (Figure  3).  The  parameterization  of  the  iono- 
grams  as  expressed  by  the  4 or  3 -bit  status  word  actually  produces 
16  (or  8)  ionograms  within  each  set  of  ionogram  data.  Printing 
the  amplitudes  without  regard  to  the  status  word  creates  a conven- 
tional ionogram.  The  electrostatic  printer  of  the  Digisonde  al- 
lows to  print  two  ionograms  side-by-side  with  128  range  bins  each. 
A set  of  80  dual  ionogram  combinations  can  be  selected  for  print- 
out. Some  interesting  examples  are: 

Positive  Doppler  - Negative  Doppler 


Low  Doppler 
Or d inary 


- High  Doppler 

- Extraordinary 


Range  Window  1 - Range  Window  2 


East  Echoes 


- West  Echoes 


Vertical  Echoes  - Non-Vertical  Echoes 


All  Echoes 


- North  Echoes 


If  the  selected  two  parameters  are  exclusive,  like  positive  and 
negative  Doppler,  the  amplitude  data  will  be  distributed  amongst 
the  ionogram  frames  according  to  the  sign  of  the  Doppler  frequency i 
blanks  will  appear  at  the  empty  FRB’s.  A special  combination  is 
the  folded  ionogram  (status  disregarded)  together  wit.',  the  status 
ionograms,  showing  the  status  values  rather  than  the  amplitudes. 

An  example  is  shown  in  Figure  4 where  the  folded  amplitude  iono- 
gram (bottom)  is  printed  together  with  the  incidence/ polarization 
status  ionogram  (top),  the  darker  numbers  indicate  the  0-component. 

For  experiments  requiring  fast  response  by  the  operator 
we  prefer  the  plasma  display  over  the  printer.  Even  though  th.e 
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electrostatic  printer  is  very  fast  (a  line  with  256  numbers  is 
printed  in  about  60  msec),  the  ionogram  does  not  become  visible 
until  after  the  paper  has  moved  out  of  the  printer.  The  plasma 
display,  on  the  other  hand,  immediately  shows  the  data  for  each 
transmitted  frequency.  The  plasma  display  and  the  hard  copy 
printout  have  cursors  inserted  indicating  the  FRB's  used  in  the 
previous  drift  measurement. 

For  M.I.'s  it  is  more  important  to  present  ionospheric 
characteristics  (Nakata  et  al,  1953)  showing  absorption  (echo  am- 
plitudes), top  E and  top  F frequencies  and  virtual  heights  as 
function  of  time.  Analog  ionosondes  have  recorded  these  charac- 
teristics on-line  for  many  years  (Bibl,  1960).  Some  Digisonde 
stations  (U.S.  Army  in  Fort  Monmouth,  New  Jersey  and  Air  Force 
Geophysics  Laboratory  in  Goose  Bay,  Labrador)  print  digital  char- 
acteristics in  real  time  (Reinisch  and  Smith,  1977).  This  fea- 
ture will  be  implemented  in  the  new  Digisonde  described  here. 

DOPPLER-DRIFT  MODE 

The  rhombic  transmitter  antenna  illuminates  a relative- 
ly large  area  in  the  ionosphere,  with  a diameter  of  several  hun- 
dred kilometers.  Echoes  return  to  the  receiver  site  from  all 
surfaces  of  adequate  electron  density  that  fulfill  the  perpendic- 
ularity condition.  A Doppler  shift  will  be  imposed  on  the  re- 
flected radio  waves  (Pfister,  1970)  when  the  ionospheric  plasma 
moves,  which  is  almost  always  the  case.  The  Doppler  shift  d is 
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determined  by  the  velocity  v of  the  reflector  and  the  radio  wave 
vector  k: 

d = (I/tt)  k*v  . 

Hence,  the  spectral  analysis  separates  different  reflection  areas 
and  it  resolves  the  angular  spread  in  each  area  (Pfister  et  al , 
1968-1975).  There  will  be  cases  where  the  same  Doppler  frequency 
will  be  generated  in  different  reflection  points;  by  using  four 
or  more  spaced  receiving  antennas  these  cases  can  generally  be 
resolved. 

Figure  5 relates  the  angular  deflection  to  the  asso- 
ciated spectral  variation.  Two  neighboring  sources  (reflection 
points)  differ  in  their  Doppler  frequency  by 
5d  = (I/tt)  k v 6(costp) 

where  i(»  is  the  angle  between  v and  k.  Using  the  angles  of  Figure 
5,  where  v was  assumed  to  be  horizontal,  we  write: 
cos^  = sinO  cos((J>-<{> ' ) 

and 


6 cosi|»  = cos9  cos((p~(P’)  66  - sin0  sin(^-^')  6^. 

When  looking  in  the  direction  of  the  velocity  vector,  i.e.  ♦ = <{>' 
the  variation  in  Doppler  becomes: 

fid  = (l/n)  k V COS0  fi0  . 

For  a given  spectral  resolution  fid  the  angular  resolution  is 
therefore: 

fi0  = If  fid/(k  V cos6)  = X fid/(2  v cos0). 

By  substituting  some  typical  values  fid  = 0.03  Hz,  X = 50  m (6 
MHz)  and  v = 100  m/s  one  obtains  for  small  angles  0 
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60  = 7.5*10"* 

corresponding  to  a source  separation  of  1.5  km  in  200  km  alti- 
tude. By  storing  the  energy  from  such  closely  spaced  sources  in 
different  spectral  channels  we  can  determine  their  location  even 
in  the  presence  of  other  strong  sources  with  different  Dopplers. 

When  looking  in  the  direction  of  v,  sources  with  60  = 0 separated 
by  a small  angle  6(|)  have  the  same  Doppler  and  can  only  be  re- 
solved by  a sufficiently  large  antenna  aperture. 

The  drift  mode  provides  the  complementary  data  to  the 
ionograms.  Up  to  six  FRB's  are  selected  for  accurate  Doppler  and 
incidence  angle  measurements.  A total  of  24  signal  channels  are 
simultaneously  processed.  Subsets  of  4 , 6,  8,  12  or  24  receiving 
antennas  can  be  selected  permitting  6,  4,  3,  2 or  1 FRB's, 
respectively. 

Receiving  Antenna  Array.  Figure  6 shows  the  antenna 
array  that  was  developed  for  DNA's  wideband  satellite  experiments 
‘ at  Kwajalein  (3°  magn.  north)  in  which  satellite,  rocket  and 

sounding  observations  are  combined  for  the  investigation  of  the 
equatorial  ionosphere.  The  aperture  of  this  25-antenna  array  is 
about  400  m and  the  minimum  spacing  is  50  m.  The  scale  size  of 
the  individual  antennas  is  2 m,  i.e.  small  compared  to  the  an- 
tenna separation.  This  is  important  in  order  to  avoid  coupling. 

The  individual  loop  antennas  are  carefully  terminated,  and  wide- 
band preamplifiers  in  each  loop  feed  the  50  Q coaxial  cables  that 
carry  the  signals  to  the  ANTENNA  SWITCH.  While  the  antennas  are 
being  scanned  the  cables  of  all  unused  antennas  are  terminated  by 
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50  in  the  ANTENNA  SWITCH;  this,  again,  is  done  to  eliminate 

coupling  by  reradiation.  Even  though  the  drift  mode  never  scans 

more  than  24  antennas  a 25th  antenna  was  added  at  the  center  of 

« 

the  array  for  the  formation  of  smaller  symmetric  sub-arrays  con- 
sisting of  only  4,  7,  or  10  antennas  (Table  5).  The  ten  inner 
antennas  actually  consist  of  crossed-loops  as  described  earlier. 
All  35  antenna  elements  (25  antennas,  ten  of  which  have  double 
polarization)  are  connected  to  the  central  ANTENNA  SWITCH  which 
either  scans  the  array  in  the  drift  mode,  or  selects  groups  of 
antennas  in  the  ionogram  mode.  Different  groups  of  antennas  are 
selected  by  simple  keyboard  instructions. 

Drift  Programs.  Observational  programs  that  are  com- 
bined with  rocket  or  satellite  measurements,  or  that  are  carried 
out  in  the  rapidly  changing  auroral  ionosphere  require  fast  deci- 
sions by  the  operator  with  regard  to  selection  of  the  "best" 
FRB's,  antenna  aperture,  spectral  resolution  and  range  and  the 
sequence  of  measurements.  Table  5 lists  the  values  of  the  opera- 
tional parameters  available  for  the  Doppler-Drift  mode.  These 
parameters  can  be  selected  via  the  front  panel  keyboard  on  the 
PROCESSING  CONTROLLER,  or  they  can  initially  be  stored  in  the 
MICROCOMPUTER  for  use  at  the  appropriate  time. 

The  maximum  number  of  data  channels  is  again  24,  as  in 
the  ionogram  mode.  Only  there  we  have  128  ranges  and  the  chan- 
nels are  used  for  selected  Doppler  lines,  polarizations,  etc.  In 
the  drift  mode  there  are,  in  general,  128  spectral  lines  and  the 
data  channels  are  divided  between  antennas  and  FRB's.  If  the 
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full  24  antenna  array  is  scanned,  only  one  FRB  is  monitored,  for 
four  antennas  six  FRB's,  etc.  Previous  drift  measurements  with 
older  Digisondes  that  were  repeated  every  five  minutes  showed 
even  in  mid-latitudes  discontinuities  between  consecutive  time 
windows  which  made  it  difficult  to  interpret  the  data.  The  new 
Digisonde  allows  for  a minimum  window  spacing  of  three  sec,  re- 
sulting, of  course,  in  a limited  spectral  resolution  of  1/2  Hz 
(Table  5). 

Discrete  Fourier  Transform.  The  receiver  output, 

IF  = 225  kHz,  is  first  logarithmically  compressed  and  then  digi- 
tized in  a quadrature  method.  The  quadrature  samples  X and  Y are 
used  directly  to  calculate  the  complex  Fourier  series  for  each 
channel  j : 

F.,  = (1/T)  I [X.(nAt)  + iY.(nAt)]  exp-i(2ii«,  nAt/T) 
jx,  ^2  2 

= (l/T)  I [X  cos6t„  * Y 
n 

* <i/T)  I [-X  sin6j„  * Yj„  cosBj„] 

where  S,  is  the  Doppler  line,  T is  the  length  of  the  time  window 
and  = 2n£  nAt/T.  While  the  initial  multiplication  of  the 

time  samples  with  the  trigonometric  functions  is  carried  out  in 
the  logarithmic  domain,  the  final  summation  is,  of  course,  per- 
formed in  the  linear  domain.  Hanning  weighting  is  applied  in  the 
frequency  domain  by  averaging  three  adjacent  spectral  lines  with 
weights  of  1-2-1.  The  result  of  the  24  channel  Fourier  analysis 
are  24  nuoiber  arrays  of  the  form 


A. 18 


1 


• • • 1 


24 


J 


where  are  the  logarithmic  spectral  amplitudes  (64  dB  range  in 
1 dB  increments,  or  32  dB  with  1/2  dB)  and  <|i.,  are  the  corre- 
spending  phases  with  a resolution  of  2n/256.  It  was  necessary  to 
develop  a method  of  calculation  that  requires  no  extra  time  be- 
yond the  last  time  sample,  since  the  time  between  windows  is  less 
than  one  second,  used  for  recording  of  the  spectra  on  magnetic 
tape  for  later  computer  processing.  Directly  following  digitiza- 
tion the  quadrature  pair  is  multiplied  by  all  128  time-adjus  ed 
values  of  cosB^^^^  and  sing^^^,  respectively,  and  the  256  spectral 
samples  thus  obtained  are  added  to  the  previous  256  samples. 

After  the  last  time  sample  the  final  four  sums  ^ X.^  cosg^^^, 

I X.  sing^^,  I Y.  cosgjj^  and  I Y.  sing,  are  suitably  combined 
n n n 

for  each  j -channel  to  form  the  real  and  imaginary  parts  of  the 


spectra . 

Incidence  Angles.  The  Doppler-drif t mode  generates  a 
sequence  of  momentary  pictures  of  the  ionosphere's  structure  and 
motions.  As  a first  step,  we  create  for  each  time  window  a sky 
map  that  shows  the  position  of  all  the  simultaneous  reflection 
points  existing  during  the  measurement,  together  with  their  in- 
tensity and  their  Doppler  frequency.  An  example  of  a sky  map  is 
shown  in  Figure  7.  We  manually  connected  the  reflection  points 
by  an  ariK)wed  line  that  pxjints  in  the  direction  of  increasing 
Doppler  frequency.  Of  course,  each  FRB  i^equires  its  own  sky  map. 
For  the  on-line  sky  maps  the  Digisondc  selects  one  FRB,  the  other 


FRB's  can  be  processed  off-line  in  a computer  using  the  tape  re- 
corded data. 


To  produce  the  angular  spectra,  the  frequency-wave- 
number-power density  is  calculated  (Sales  et  al,  1975) 


£k 


: D 


where  F..  is  the  spectrum  for  antenna  j , k is  the  wave  vector  and 
J *’ 

specifies  the  looking  angle,  and  a.  is  the  antenna  distance  vec- 
tor; is  a weighting  factor  (Varad  and  McComish,  1975)  which 

optimizes  the  array  performance.  For  the  on-line  generation  of 
the  sky  maps  in  the  Digisonde,  the  calculations  use  four  anten- 
nas, and  since  in  that  case  antenna  weighting  does  not  offer  any 
advantage,  we  set  =1.  The  power  density  can  then  be  re- 

written as 
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hk  = h Ajt  «=kp  i (♦jj-k-tj)!’ 

where  we  substituted  F . . = A . „ exp  i<j) . - . For  each  looking  angle, 
specified  by  k,  the  spectrum  must  be  calculated.  The  sky  map 

is  composed  of  40  40  points,  which  means  that  1600  spectra  must 

be  calculated  in  a tew  seconds.  Even  a large  computer  would  be 
hard  pressed  to  master  this  task.  We  developed  special  micro- 
computer controlled  hardware  that  calculates  the  1600  power  spec- 
tra in  about  three  seconds.  It  is  practically  impossible  to 

print  out  the  complete  spectrum  for  each  individual  k.  Instead, 
we  determine  for  each  Doppler  component  its  position  in  the  sky 
(indicated  by  maximum  power).  If  more  than  one  Doppler  component 
occupies  the  same  point,  the  one  with  the  largest  amplitude  is 
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selected.  Both  the  power  and  the  number  of  the  Doppler  line  are 
presented. 


1 
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Ionospheric  Modelling.  Once  the  ionospheric  reflection 
points  are  found  one  can  proceed  to  fit  ionospheric  models  to 
these  points.  The  mathematics  and  the  computer  software  for  this 
problem  are  presently  being  developed  on  the  basis  of  assuming  a 
two  dimensional  ionospheric  density  structure  described  by  two 
wave  vectors  and  ^2  that  lie  in  a tilted  plane: 

z = Zq  + ax  + by  + A^costflj^*  (r-r^)  ] + A2Cos[ft2  * ^£“£2  ^ ^ 

The  constraint  in  the  fitting  procedure  is  the  condition  of  per- 
pendicularity at  the  observed  reflection  points.  Preliminary  re- 
sults for  the  location  of  the  reflection  points  by  a one-dimen- 
sional structure  were  derived  earlier  (Bibl,  1974),  but  our  ob- 
servations seem  to  suggest  isodensity  surfaces  that  are  modulated 
in  two  dimensions. 

The  method  of  calculating  the  three-dimensional  drift 
vector  from  the  distribution  of  Doppler  frequencies  d = d(k)  was 
described  in  an  earlier  paper  (Bibl  et  al,  1975).  The  calculated 
drift  vectors  together  with  the  sequence  of  sky  maps  will  reveal 
the  small  and  large  scale  motions. 

CONCLUSION 

f 

The  Digisonde  128PS  represents  a systematic  approach  to 
the  measurement  of  the  ionosphere  with  radio  waves.  All  wave 
signatures  are  exploited  to  retrieve  information  on  the  structure 
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and  the  dynamics  of  the  reflecting  medium.  The  large  variety  of 
observational  programs  is  conveniently  organized  in  a set  of  sub- 
routines that  make  fast  program  changes  possible.  Selection  of 
individual  operational  parameters  is  internally  controlled  and  a 
default  algorithm  corrects  the  selection  if  an  illegal  parameter 
combination  was  requested.  The  phase-coherent  spectral  integra- 
tion provides  good  signal-to-noise  ratio  even  in  environment;; 
with  high  interference,  or  in  cases  where  a small  transmittiig 
antenna  or  limited  power  is  required.  An  automatic  f requeue/  se- 
lection scheme  tests,  in  the  ionogram  mode,  three  closely  spaced 
frequencies  around  the  nominal  frequency  for  the  lowest  interfer- 
ence before  the  start  of  integration.  Digital  automatic  gain 
control  further  improves  the  appearance  of  the  ionograms.  A time 
dependent  gain  control,  programmed  for  the  entire  year,  adjusts 
the  receiver  gain  for  day,  twilight  and  night  condition^  in  6 dB 
increments . 

We  believe  that  this  instrument  will  create  new  knowl- 
edge for  the  physical  explanation  of  spread  F conditions,  travel- 
ling ionospheric  disturbances,  the  high  latitude  trough,  scintil- 
lations, as  well  as  radio  communication  problems  involving  the 
ionosphere. 
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APPENDIX  B 

IONOSPHERIC  MONTHLY  MEDIAN  CURVES 
FOR  GOOSE  BAY,  MARCH  1975  - FEBRUARY  1976 
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4 10  16  22  UT  04 

LIMITING  FREQUENCY  = 3MH2  

LIMITING  FREQUENCY  =5MH2 

LIMITING  FREQUENCY  =7MHi  

f 

i 

i 

i 

f 


B.2 


APRIL  1975 
GOOSE  BAY,  LABRADOR 


LIMITING  FREQUENCY  = 3MHi  

LIMITING  FREQUENCY  =5MHi 

LIMITING  FREQUENCY  =7MHi  


MAY  1975 

GOOSE  BAY,  LABRADOR 


LIMITING  FREQUENCY  =3MHz 
LIMITING  FREQUENCY  =5MHz 
LIMITING  FREQUENCY  =7MH« 


PERCENTAGE  OF  TOTAL  TIME  DURING  WHICH 
toEj  IS  GREATER  THAN  THE  LIMITING  FREQUENCY 

FIGURE  6 
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JUNE  1975 


GOOSE  BAY,  LABRADOR 


4 10  16  22  UT  04 

LIMITING  FREQUENCY  = 3MH2  

LIMITING  FREQUENCY  =5 MHi 

LIMITING  FREQUENCY  =7MHi  


PERCENTAGE  OF  TOTAL  TIME  DURING  WHICH 
foE,  IS  GREATER  THAN  THE  LIMITING  FREQUENCY 


FIGURE  8 
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JULY  1975 

GOOSE  BAY,  LABRADOR 


LIMITING  FREQUENCY  = 3MHi  

LIMITING  FREQUENCY* 5 MHx 

LIMITING  FREQUENCY  =7MH*  


PERCENTAGE  OF  TOTAL  TIME  DURING  WHICH 
foE,  IS  GREATER  THAN  THE  LIMITING  FREQUENCY 

FIGURE  10 
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SEPTEMBER  1975 
GOOSE  BAY,  LABRADOR 


LIMITING  FREQUENCY  = 3MHz  

LIMITING  FREQUENCY  =5MH2 

LIMITING  FREQUENCY  =7MHi  


PERCENTAGE  OF  TOTAL  TIME  DURING  WHICH 
foE,  IS  GREATER  THAN  THE  LIMITING  FREQUENCY 


FIGURE  14 


B.  14 


B.15 


r 


OCTOBER  1975 
GOOSE  BAY,  LABRADOR 


LIMITING  FREQUENCY  = 3MHz 
LIMITING  FREQUENCY  =5MHi 
LIMITING  FREQUENCY  =7MHi 


PERCENTAGE  OF  TOTAL  TIME  DURING  WHICH 
toE,  IS  GREATER  THAN  THE  LIMITING  FREQUENCY 

FIGURE  16 
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NOVEMBER  1975 
GOOSE  BAY,  LABRADOR 


LIMITING  FREQUENCY  = 3MHz 
LIMITING  FREQUENCY  =5MH2 
LIMITING  FREQUENCY  =7MHz 


PERCENTAGE  OF  TOTAL  TIME  DURING  WHICH 
foEj  IS  GREATER  THAN  THE  LIMITING  FREQUENCY 


FIGURE  16 
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FIGURE  22 
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I.  INTRODUCTION 


This  manual  is  primarily  a users'  guide  to  digital 
ionograms.  It  is  not  meant  to  be  sufficient  unto  itself;  it 
must  be  read  with  the  U.R.S.I.  Handbook  of  lonogram  Interpre- 
' tation  and  Reduction.  Users  of  evaluated  data  as  well  as  raw 

and  processed  digital  ionograms  may  find  additional  informa- 
tion furnished  in  this  guide. 
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II.  IDENTIFICATION  OF  TRACES 


I 

i 

I 

I 


f 


We  are  concerned 

E Region 
F Region 


with  two  regions: 

90-150  km 
150-750  km 


E-Region 

Normal  E - A normal  layer  whose  maximum  frequency 
increases  with  cos  x (solar  zenith  angle).  This  layer  is 
thick  and  shows  significant  retardation  (cusp)  near  the 
layer's  maximum  ionization  density  (Figure  1). 


Sporadic  E - An  abnormal,  and  mostly  thin,  layer 
caused  by  localized  ionization;  this  layer  can  have  a very 
short  lifetime.  The  Es  layer  is  classified  into  many  types 
dependent  on  its  shape  and  locality. 


} 

i 


Auroral  - diffuse;  usually 
nonblanketing;  flat  or  slowly 
increasing  height  with  fre- 
quency with  stratified  traces 
in  it  which  vary  rapidly  with 
time;  most  common  Es  type  in 
Goose  Bay.  (See  also  URSI 
Handbook,  Fig.  4.29.) 

Cusp  - Es-ionization  located 
near  the  height  of  maximum 
E layer  ionization;  usually 
decreasing  delay  with  in- 
creasing frequency. 
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Flat  - classification  used 
only  at  night,  when  a thick 
E-layer  is  not  usually  ob- 
servable. At  other  times 
this  Es  will  be  classified 
as  low. 

High  - ionization  located 
above  maximum  height  of  E. 

Low  - flat  Es-,  ionization 
usually  located  below  mini- 
mum height  of  E-,  similar  in 
appearance  to  F type  Es. 

e2,uatorial  - similar  in 
appearance  to  Auroral  Es,  yet 
caused  by  different  processes. 

Retardation  - thick  Es  layer 
similar  in  appearance  to  E. 
Retardation  in  high  frequency 
end . 

Slant  - similar  to  Auroral, 
the  S type  usually  develops 
from  another  Es. 


The  most  often  used  type  classifications  for  Goose 
Bay  are  A,  C,  F and  L.  The  H and  S types  may  be  labeled  as  A 
because  they  often  have  the  same  characteristics. 
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Night  E - A thick  layer  which  shows  the  same  I'etar- 
dation  trends  as  a normal  E yet  has  no  r'elation  to  x-  The 
night  E is  very  special  in  that  it  is  recorded  as  both  an  E 
and  Es  simultaneously.  Night  E is  usually  found  to  resemble 
one  of  the  four  configurations  described  in  Figures  2 through 
I 5.  Note:  Frequently  foE  is  more  correctly  determined  by 

' fminF  (minimum  frequency  of  F layer)  than  by  the  inside  cusp 

of  the  E trace. 

F-Region 

Normal  F - Some  F region  ionization  Ls  always  pre- 
sent. During  the  daytime  (when  x affect  the  ionization 

structure)  we  have  a separation  into  two  layers,  FI  and  F2 
(shown  in  Figure  1).  The  FI  layer  shows  a frequency  pre- 
dictability similar  to  the  normal  E.  The  maximum  frequency 
of  the  FI  foliC'WS  the  relation 

foFl  E foE  X 1.5  . 

Transitoria  or  Intermediate  Layers  - Fast  moving 
stratifications  which  are  embedded  in  the  F ionization  or 
exist  as  a layer  between  the  E and  F region.  Intermediate 
layers  are  often  referred  to  as  E2  layers. 

Oblique  F (Figure  6)  - In  regards  to  vertical  in- 
cidence data  U.R.S.I.  neglects  all  oblique  readings;  however, 
f 3r  our  arctic  studies  these  values  are  important.  Oblique  F 
echoes  always  have  a longer  delay  than  normal  F;  during  a cer- 
tain time  period  the  oblique  trace  can  move  in  and  become  an 
overhead  echo.  Since  oblique  traces  are  often  only  fragments, 
it  is  not  always  possible  to  differentiate  between  oblique  F 
and  long  range  oblique  Es. 
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h'E=  119  K ; h'Es=  119  K ; h'r=  272EA 
hpF2=  R ; TYPE  Es=  R (km) 

fmin=  1.3  ; foE=  4.5  K ; fbEs=  4.5  ; 
foEs=  4.5  K ; fxE=  fxLs=  5.4  K ; 
foF2=  5.4UR  (MHz) 
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h'E=  119  K ; h'Es=  119  K ; h'F=  308  ; 
hpF2=  F ; TYPE  Es=  R (km) 

fmin=  0.7  ; foE=  1.5  K ; fbEs=  1.5  ; 
foEs=  1.8  K ; for2=  A (MHr.) 


h'E=  114  K ; h'Es=  114  K ; h ' F=  A;  hpF2=  A 
TYPE  Es=  R 

fmin=  1.6  ; foE=  5.2  K ; foEs=  5.9  K ; foF2= 
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III.  TAKING  OF  READINGS 


Each  trace  has  up  to  three  magneto-ionic  components 
(Figure  7)  which  may  be  visible  simultaneously:  the  ordinary 
(0),  extraordinary  (X),  and  Z mode  (Z).  The  frequency  sepa- 
ration between  Z and  X is  equal  to  fH  (gyrof requency ) . For 
Goose  Bay  the  gyrofrequency  is  approximately  1.52  MHz  at 
100  km.  Therefore  the  equation  for  frequency  mode  is 

fo  - fz  = 1/2  fH 
fx  - fo  = 1/2  fH 

thus , 

fx  - fz  = fH 

The  normal  E and  F layers  are  considered  monotonic; 
this  means  that  the  curve  of  ionization  density  vs  range  is  a 
continuous  trace,  therefore  foE  <_  fminF ; foFl  ^ fminF2. 

For  Digisonde  ionograms  there  are  three  criterion 
for  a trace: 

1)  the  amplitude  must  be  at  least  6 dB  greater  than 
background  noise  (Note:  be  very  sure  what  you  see 
is  not  merely  a spike  on  a noisy  frequency  line); 

2)  the  trace  has  a finite  thickness  - in  K = 3 (verti- 
cal ionograms)  the  pulse  width  is  4 height  bins,  in 
K = 7 (backscatter ) the  pulse  width  is  2 height 
bins ; 

3)  the  virtual  height  of  the  echo  is  the  highest  range 
bin  where  the  amplitude  is  20  dB  (five  units)  below 
the  average  maximum  amplitude. 
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The  most  significant  readings  are: 

h'  - minimum  virtual  height 

h'E;  h'Es;  h'F;  h'F2 

When  there  is  only  one  F layer  its  height  is  taken 
as  h'F. 

fo  - maximum  frequency  of  ordinary  trace 

foE;  foEs;  foFl;  foF2 

When  there  is  only  one  F layer,  its  maximum  critical 
frequency  is  foF2  except  in  G-effect  cases  when  the 
ordinary  trace  is  not  easily  discernible  as  in  the 
case  of  spread  E or  spread  F conditions;  foE  or  foF2 
has  to  be  read  from  the  inner  edge  of  the  spread  E 
or  spread  F pattern.  (See  URSI  Handbook,  page  41, 
and  Fig.  2.11,  and  Fig.  12,  this  report.  For  deter- 
mination of  M(3000)  under  spread  F conditions  see 
URSI  Handbook  Fig.  3.14.) 

ftEs  - top  frequency  of  Es  - whfn  0 and  X modes  are  not 
distinguishable. 

fmin  - the  minimum  frequency  at  which  we  first  see  contin- 
uous echoes. 

fbEs  - the  frequency  at  which  the  Es  becomes  transparent  or 
nonblanketing . 

hpF2  - parabolic  height  F2.  The  height  of  the  F2  layer 
when  the  frequency  is  .834  x foF2. 

frEs  - retardation  frequency  of  Es,  as  shown  in  Figure  11. 


Type  Es  - classification  of  Es  according  to  its  appearance. 
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IV.  QUALIFYING  OF  READINGS 


If  there  is  any  doubt  or  further  explanation  is 
needed  to  clarify  a reading,  qualifying  and  descriptive  let- 
ters are  added  to  the  numerical  value.  Qualifying  letters 
cannot  exist  without  descriptive  letters,  yet  the  descriptive 
letters  may  be  used  alone.  In  the  case  of  a probable  error 
> 15%  just  a descriptive  letter  is  used  with  numerical 
reading . 


Qualifying  Letters 

With  respect  to  the  actual  reading  the  real  value 

would  be: 

D - greater  than  5%  < error  < 10% 

E - less  than  5%  < erroi"  < 10% 

I - interpolated 

J - derived  from  X component 

U - approximately  2%  < error  < 5% 

If  the  error  is  < 2%  no  qualifier  is  necessary. 


Descriptive  Letters 

Measurements  influenced  or  impossible  because  of : 

A - blanketing  (Figures  8 to  10) 

B - absorption  (Figure  10) 

C - technical  difficulties 
D - upper  frequency  limit 
E - lower  frequency  limit 
F - spread  (Figure  8,  Figures  11  and  12) 

G - low  ionization  density  (Figure  10) 

H - stratification 

K - night  E (Figures  2 to  5 , Figure  12) 
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L - no  cusp  (Figure  11) 

M - 0 and  X not  distinguishable 
0 - ordinary  component 
Q - range  spread  (usually  use  F) 

R - attenuation  near  maximum  frequency  (Figure  13) 
S - radio  interference 
T - derived  from  other  readings 

V - forked  trace 

W - outside  height  range 
X - derived  from  X trace 

Y - F lacuna 

Z - derived  from  Z trace 


fmin 


fmin 


fmin=  1.1  ; foF2=  6.1 


fmin=  1 . 3 


FIGURE  13 


V.  SPECIAL  EXAMPLES 


VI.  PROCESSING 


There  are  several  modes  of  processing  which  have 
f been  used  to  present  the  Digisonde  ionograms  in  diurnal 

displays . 


A.  Compressed  Data  (Figure  14) 

h'(t)  - each  height  bin  was  summed  along  all  fre- 
quencies of  one  ionogram  and  the  array  of  128  summed  height 
bins  was  then  normalized  to  63  using  the  relation 


Anew  = ^ ) 

max  min 


(A. -A  . ) 
1 min 


ftF(t)  - maximum  amplitude  of  height  bins  > 156  km 
for  each  frequency  was  saved  and  a display  of  frequency  vs 
amplitude  was  made. 


ftE(t)  - maximum  amplitude  of  height  bins  < 156  km 
is  saved  and  displayed. 


B.  Reconstituted  Ionograms  (Figure  15) 

Using  a process  of  matching  a pulse  filter  along 
each  height  bin  of  a certain  frequency,  six  significant  traces 
are  extracted  from  each  ionogram.  These  traces  may  then  be 
displayed  as  the  reconstituted  ionogram  which  has  the  informa- 
tion of  amplitude,  height  and  spread  of  the  original  ionogram. 
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The  method  which  can  create  the  reconstituted  iono- 
grams also  produces  the  many  displays  shown  in  Figures  16  and 
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FILTERED  VERTICAL  lONOGRAM 
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